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ABSTRACT 



We perform three-dimensional gas-dynamical simulations and show that the 
asymmetric morphology of the blue and red-shifted components of the outflow at 
hundreds of astronomical units (AU) from the massive binary system r] Carinae 
can be accounted for from the collision of the free primary stellar wind with the 
slowly expanding dense equatorial gas. Owing to the very complicated struc- 
ture of the century-old equatorial ejecta, that is not fully spatially resolved by 
observations, we limit ourselves to modelling the equatorial dense gas by one or 
two dense spherical clouds. Because of that we reproduce the general qualitative 
properties of the velocity maps, but not the fine details. The fine details of the 
velocity maps can be matched by simply structuring the dense ejecta in an appro- 
priate way. The blue and red-shifted components are formed in the post-shock 
flow of the primary wind, on the two sides of the equatorial plane, respectively. 
The fast wind from the secondary star plays no role in our model, as for most of 
the orbital period in our model the primary star is closer to us. The dense clouds 
are observed to be closer to us than the binary system is, and so in our model 
the primary star faces the dense equatorial ejecta for the majority of the orbital 
period. 



Subject headings: stars: mass loss — stars: winds, outflows 
other — stars: individual (Eta Car) 



stars: variables: 



INTRODUCTION 



A picture is emerging in recent years according to which binary interaction plays a major 



role in the evolut ion of very massive stars. The main 



processes involved are mass transfer 



and merger (e.g., iKashi fc Sokeiil2010al : ISoker fc Kashill2012al : ISana et al.l 120121 ). The best 
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example is rj Carinae, one of the most massive binary systems in our galaxy, composed of a 
primary Luminous Blue Variable (LBV), and a somewhat evolved m assive 0-star secondary . 
The estimated masses are Mi ~ 120 - 2OOM0 and M2 ~ 60-90M^ kashi fc Sokeilbninah . 



respectively. The orbital period is five and a half years (iDamineli et al.ll2008l ). and the orbit 
is highly eccentric, e ~ 0.9. For several weeks every periastron passage the binary system 
experiences a strong interaction, possibly with some mass transfer. Along the rest of the 
orbit the winds from the two stars collide with each other. The colliding winds determine 



comes 


directly from the post-shock regions 


f t 


he two winds (Corcoran 2005: 


Corcoran et al. 


2010 




Parkin et al. 


201 1|; 


Akashi et al. 2006 




2012; 


Moffat & Corcoran 


200c 


); 


Henlev et al. 


2008 


Okazaki et al. 


2008; 


Pittard & Corcoran 


2002; 


Pittard et al. 1998 


Behar et al. 2007). 



The binary interaction was much stronger during the two eruptions in the nineteenth 
century, when the primary is thought to have been in an unstable state. Huge amounts of 
mass were lost during the 1837.9 - 1858 Great Eruption (GE; see iDavidson fc Humphreys 
19971 for a review on the GE and some other observed properties of t] Car). The bipolar 
Homunculus nebula surrounding the system an d containing most t he e lected mass is the 



most prominent structural feature from the GE ( Gomez et al 



2006 



2006 



2010 



Smith &: Owocki 



Smith et al.l 120031 : ISmithfc Ferlandl 120071 ) . Some non-negligible fraction of the mass 



l ost by the primary s tar was accreted by the secondary star during these twenty years 
( iKashi fc Sokerll2010al ). The accretion rate was higher near periastron passages, accounting 
for observed luminosity peaks during the GE. 

Both the GE and the weaker 1887.3 - 1895.3 Lesser Eruption (LE) that followed it are 
the sourc e of most of the circurn binary medium (CBM). The LE was a much less en ergetic 
erupt ion f lHumphrevs et al.lll999l ) and only 0.1 - IM0 were ejected from the primary ([Smith 
2OO5I ). Material ejected during the LE is thought to be the source of the de nse inhomo- 



geneo us blobs and the gas in their surroundings closer to th e binary system (Smith 



et al. 



1986 



2004). These blobs are k nown as the Weigelt blobs (WBs; Weigelt &: Ebersbergerl 
Hofmann fc Weigeltl Il988l ) . The WBs are concentrated on the NW side of r] Car that is 



closer to us, but there is also dense mater i al in other direc t ions fe.g.. ISmith fc Gehrzl 12000 



Dorland et al.ll2004t IChesneau et al.l 120051 IGuII et al.l 120091 : iMehner et al.l 120101 ). The inter- 
action of the primary wind with the WBs environment (WBE) is the subject of our study. 

From HST/ST IS spectra of Fe II, [Fe HI], [Ni II] and [Ni HI] lines from WBs C and D, 
Smith et al.l (120041 ) found the radial velocity of the WBs to be ^ 40 km s~\ and concluded 
that the WBs were originated in the LE. On the other hand, iDorland et al.l (120041 ) presented 
HST observations of WBs C and D together with simulations of their propagation and 
concluded that they formed during the period 1910 - 1942, preferably during a brightening 
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event which took place in 1941. For the purpose of this paper it is immaterial when the 
dense WBE was ejected, but the uncertainty should be kept in mind. 

The most controversial property of the binary system is the direction of the semi-major 
axis in the equatorial plane, the so called longitude angle u (the argument of periapsis). 
It is defined such that u = 270° for the case where the secondary is closest to us at apas- 
tron passages, and it is furthest away from us during periastron passages. For u = 90° the 
secondary is furthest from us at apastron passages, while it is closest to us at periastron 



a literally opposite view, with one camp arguing for u ~ 270° (e.g., 


Moffat & Corcoran 




2009; 


Okazaki et al. 


2008; 


Henlev et al. 


2008; Hamaeuchi et al. 


2007: 


Nielsen et al. 2007 


Gull et al. 2011: iMadura et al. 2011. 201 


2: Madura & Groh 2012: 


Parkin et al. 2009). while 


the other camp arguing for oj ~ 90° (e.g.. 


Abraham & Falceta-Goncalvesll2007: Falceta-Goncalves et al. 


2005; 


Kashi & Soker 


2008, 


2009c). 



The quest for the longitude angle is not a small technical curiosity. Its value is connected 
to the type of interaction between the winds of the two stars and between the winds and the 
CBM. For example, if the secondary is towards us near periastron passages, u ~ 90°, then 
the several- weeks decline in the X-ray emission occurring around periastron passages must 
be acc ou nted for by accreti on of the primary wind by the secondary star (IKashi fc Soker 
2009al lbl; | Akashi et al.ll2012l ). This is because the dense primary wind cannot account for 
much X-ray obscuration. This in turn, implies that accretion was much more vigorous 
during the GE, with implications to the energy source of the GE. Namely, a value of u ^ 
90° su pports the claim for accretion as the extra energy source of the GE (IKashi fc Soker 
2010ar). hinting to an accre tion event as the energy source for other transient outbursts 
kashi fc Sok3l2010bl . Eoill ). 



In this paper we examine the interaction of the primary wind with a simple model of 
the WBE. In section [2] we describe the flow structure and previous works. The numerical 
setup is described in section [31 The interaction of the primary wind with the WBE model 
is described in section IH and a comparison between our simulation results and observations 
is conducted in section [51 Our short summary is in section [61 



THE VELOCITY STRUCTURE AT HUNDREDS OF AU 



The recent dispute over the value of the longitude angle is centered ar ound the in t erpre- 
tation of the velocity structure at hundreds of AU from the binary system. iGuU et al.l (l2009l ) 
presented observations of velocity structures of broad (~ 500 km s~^) [Fe II] emission line 
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which extend to ~ 1600 AU from the binary system, and [Fe III], [Ar III], [Ne III] and [S III] 
l ines w hich extend up to a projected distance of only ~ 700 AU from NE to SW. iGuU et al. 



(120091 ) observed all those forbidden lines disappearing during the spectroscopic event - the 
several weeks period near periastron passages when large variations in emission and absorp- 
tion properties occur. iMadura et al.l (120121 ) took HST/STIS spectra of [Fe III] emission lines 
from slits in different position angles through the central sour ce of r? Car, ^yhile t he low 
velocity component of these lines was associated with the WBs, iMadura et al.l ( 120121 ) inter- 
preted the high velocity components of these lines a s bein g formed by winds collision gas as 
it flows away from the binary system. iGuU et al.l ( 120111 ) presented more HST/STIS high- 



ionization forbidden-line observations taken after the 2009 spectroscopic event, this time in 
the form of complete intensity maps. 

Based on comparison with SPH simulations and under the assumption (which we dis- 
pute in the p resent paper) tha t the lines orig i nate b y the colliding - winds ma terial as it 



flows outward, IGuU et al.l (]2009[ ) , Madura et al.l (120111 ) , Madura et al.l tom and IGuU et al. 



(120111 ). concluded that the high velocity components of the forbidden- line observations fit 
an orientation in which the secondary star is in the direction of the observer for most of the 
binary orbit. Namely, the longitude angle (the argument of periapsis) is a; ^ 270° according 
to their interpretation. 



Mehner et al.l (120101 ) also presented HST/STIS spectra of high ionization lines ([Fe III] 
and [N III]) around rj Car extending over more than a full cycle. They also found that most of 
the emission arrives from the WBs region. Moreover, they found a variation of the intensity 
of the lines over the entire orbital cycle, in addition to the fast incre ase before per i astron 
that was followed by the typical decrease of the spectroscopic event. Mehner et al.l (120101 ) 
identified a fast blue-shifted component of the high ionization lines that appear s concentrated 



near t he stars and elongated perpendicular to the bipolar axis of the system. Mehner et al. 



(j2010l ) suggested that this component is related to the equatorial outflow and/or to dense 
material known to exist along our line of sight to the system. 



Madura et al.l (120121 ) performed 3D numerical simulations of the colliding stellar winds, 
but completely ignored the slow dense eiec ta in the equatori al pl ane - the WBE. The high 



resolution observations reported recently by lCuU et al.l ( 1201 ll ) and Madura et al.l (2012) that 



extend to ^ 0.^^5 from the binar y system overlap with the WBE as mapped by 



Chesneau et al. 



mm and 



Madura et al. 



( 2005 ). ISoker fc Kashil (1201 2br) r eexamined the observations and interpretation of iGuU et al. 



(l2012h 



and concluded that the interaction with the WBE cannot 



does not hold. Instead, 



be ignored when n ao deling emissi o n from that region. Hence, ISoker fc Kashil (l2012bl ) claimed, 
the conclusion of Madura et al.l (12012 ) that u ~ 90° cannot account for the velocity map 



Soker &: Kashil (l2012b| ) suggested that the contribution of the winds 
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collision gas to the observed velocity structure at hundreds of AU is small. The main 
contribution to the asymmetric flow structure is of the freely expanding primary wind as 
it collides with the WBE. Because the WEE is closer to us, and for most of the time the 
secondary is away from us for u ^ 90°, the secondary wind does not intercept the primary 
wind flowing towards the WBE. 



In the scenario proposed by ISoker &: Kashil f l2012bh the primary stellar wind collides 



with the slow dense material in the WBE and goes through a shock wave. The post-shocked 
gas reaches a temperature of ~ 10^ K but rapidly cools to ~ 10^ K and is compressed. This 
warm dense gas is expected to become a source of some emission lines, including forbidden 
lines. Soker Sz Kashi further showed that the intensity and location of the different 



velocity components of the [Fe III] A4659 extend ed forbidden line can b e explained by the 



interaction of the primary wind with the WBE. ISoker fc Kashil ( l2012bh made a schematic 



drawing of the flow structure but did not conduct any numerical simulations of the flow. 



In Figured] we present the observed inner region taken from Figure 5 of IChesneau et al. 



( 120051 ). It is evident that the region surrounding the dense WBs, the WBE, around hundreds 
of AU on the NW direction (upper right in the image) is quite complicated. The exact 
structure is not known as we see material in projection. It is assumed that the dense gas is 
near the equatorial plane that is inclined to line of sight by 42°. In the present paper we limit 
ourselves to show that the general velocity structure can be explained by the interaction of 
the primary wind with the WBE, as we cannot reproduce accurately the exact density near 
the equatorial plane. It will be clear that by rebuilding the WBE in the numerical code we 
can reproduce the flow structure. A big uncertainty is the structure of the dense gas and its 
velocity when it was formed a century ago. Another uncertainty in reproducing the image is 
the emissivity at each point due to the lack of knowledge of the ionizing radiation flux from 
the secondary. Indeed, the velocity maps change over the orbital period. 

Considering the above, we take very simple WBE structures, either one or two spherical 
dense blobs (clouds). More details are given in section El 



NUMERICAL SETUP 



The simulations a re performed using the high-resolution multidimensional hydrodynam- 
ics code FLASH 4.0b fiFrvxell et al.ll2000l ). We employ a full 3D uniform grid (all cells have 
the same size) with Cartesian {x,y, z) geometry. We define the numerical x axis to lie along 
the line of sight. 



The numerical y and z axes are defined so that the center of the primary star, that is 
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Fig. 1. — The dense Weigelt bl obs (WB) and th e ir env ironment, the WBE, within hundreds 
of AU from r] Car, taken from IChesneau et aLl (120051 ). The image is a deconvolved image 
of the emission line Pfund 7. At the distance of 77 Car each arcsec is 2300 AU. WBs B, C 
and D reside in the equatorial plane, inclined to the line of sight by 42°. In the inset we 
show the density in the equatorial plane of one of our simulations (see Section H]). The inset 
corresponds to the square marked around the center of the image. 



the source of the wind, the observer and the centers of the clouds that form the WBE in our 
model are all in the x — y plane {z = 0). The WBE model is symmetric about the z = 
plane. The z axis and the vector from the primary star to the center of the large WBE cloud 
define a plane that corresponds to the equatorial plane of the binary system. The grid size 
is 3000 AU X 3000 AU x 3000 AU in the x,y,z directions, with 128 equal-size cubical cells 
along each axis (single cell size is ~ 24 AU). The wind of the r] Car primary star is modelled 
by a radial outflow from the central region of the grid (the inner 100 AU), with a mass loss 
rate of Mi = 3 x 10~^ Mq yr~^ and velocity of vi = 500 km s~^. 

The WBE has been ejected by the system more than a century ago. However, some 
dense equatorial outflow might have been continued for tens of years. It is impossible to 
know the exact mass loss history and mass loss rate that led to the formation of the WBE. 
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Considering that our goals are (a) to show that the interaction of the primary wind with the 
WBE cannot be neglected, and (b) to show that this interaction can in principle give the 
velocity maps at hundreds of AU from the system, we consider two very simple models. In 
those models the WBE is composed either from one or two spherical clouds that start with 
zero velocity. A more realistic value would be to start the WBE with some radial outflow 
velocity of ~ 40 km s~^. Another uncertainty is the history of the primary wind, that could 
have been weaker until recently, or could have been stronger in the equatorial plane until 
recently. It is hard to know. We simply base our two models on the present existence of 
the WBE and its location. Considering these, we run our simulations for ~ 50 yr, until the 
wind launched from the center of the grid fills the entire grid. That our WBE started with 
zero velocity rather than 40 km s~^ should be kept in mind while comparing our results with 
observations. 

The one-cloud structure is composed of a dense spherical cloud centered 700 AU from 
the center of the primary, as depicted in the inset of Figure [T] and in Figure [H The mass 
of the cloud is Mwbei = 0.2Mq, and the number density in the center of the cloud is 
nwBE = 10^ cm~^ (pwBE = 1.03 x 10~^^ g cm~^). We assume that the cloud extends to 
a radius of 600 AU, and its density profile is given by pcioudi = Pwbe x ( 70 au+twhei '^ ^'^^ 
< fwBEi < 600 AU, so that the initial density at twhei = 600 AU is ~ O.lpwBE- 

In the two-clouds structure a smaller cloud with a radius of rwBE2 = 200 AU is added 
to the one-cloud structure in the x — y plane, at a distance of 300 AU from the primary 
star, and at an angle of 30° from the equatorial plane (or 72° from the line of sight, which is 
the positive x axis in our simulation). The mass of the smaller cloud is Mwbe2 = O.O6M0, 
but we take it to be denser than the large cloud. The density in the center of the smaller 
cloud is 2 X Pwbe, and its density profile is given by pcioud2 = 2 x pwBE x ( 200 au4-wbe2 ) 
for < rwBE2 < 200 AU, so that the initial density at rwBE2 = 200 AU is ^ Pwbe- The 
location and properties of the second cloud are chosen almost arbitrarily. The sole goal is 
to demonstrate that small changes in the WBE structure can have noticeable effect on the 
velocity maps. As we do not know the exact structure of the WBE, at this stage we limit 
ourselves to qualitative study. 

The initial temperature of the primary wind and the WBE is Tj = lOOOK. We add 
ra diative cooling for T > lO^K , using the solar composition cooling function A(r), as given 



by [Sutherland fc Dopital (119931 ). We impose outflow boundary conditions on all 6 sides of 



the simulation box. The initial numerical setup is drawn schematically in Figure [2l 
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Fig. 2. — The initial setup of the numerical grid for the one-cloud model. A 2D plane of the 
3D computational domain is shown, where the third axis is —1500 AU < z < 1500 AU. The 
location of the second cloud in the two-clouds model is marked by the green dashed circle. 
The equatorial plane of the binary system, which is taken to be where most of the WBE 
mass is, is defined by the z axis and the line from the the center of the grid to the center of 
the large WBE cloud. The intersection of this plane and the x — y plane is marked by the 
white dotted line. The observer is in the positive direction of the x axis in the z = plane. 



4. WIND- WBE INTERACTION 

We present our results after 50 years of simulation, when the post-shock primary wind 
reached the boundary of the simulation grid. In Figure [3] we present the density and velocity 
maps in the x — y plane of the one-cloud model for the WBE (see Section [3] and the inset 
of Figure [1]). Regions that contribute to the red, blue, and zero Doppler shifts are marked. 
The primary wind hits the WBE and a shock is formed. At the apex, where the flow hits the 
shock wave at a straight angle, the immediately post-shock temperature is calculated to be 
3 X 10® K. The immediate post-shock temperature decreases away from the apex because the 
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Density 
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Fig. 3. — The simulated interaction between the r] Car primary wind and the WBE in the 
X — y plane. The density scale is logarithmic in units of g cm~^. The observer is on the 
right side of the figure looking left. The WBE is closer to the observer than the central 
binary system is, and inclined to the line of sight by 42°. The post-shock material moving 
towards the observer in the dense- warm region is encircled by a blue line, while the post- 
shock material moving away from the observed is encircled by a red line. The WBE itself 
constitutes most of the zero- velocity component. The velocity vectors of the post-shock 
material and the initial wind are marked by arrows (the different colors of the arrows are for 
clarity of presentation only). Velocity arrows of the pre-shock primary wind are not plotted, 
as the flow is a simple radial outflow. The length of the arrows is proportional to the velocity, 
with a maximum at 500 km s^^. 
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shock becomes more oblique. In any case, the post-shock gas very rapidly cools radiatively 
to a temperature of T ~ 10^ K, and a temperature of ~ 10® is not seen in the temperature 
map. Actually, the post-shock gas moves a distance of less than a cell size before it cools to 
T ~ 10^ K. A warm and dense gas now occupies the post-shock regions. It is this gas that 
is the main source of the radiation observed in the velocity maps according to our model. 
As can be seen in Figure [31 the freely flowing primary wind is diverted from its initial radial 
direction as it passes through the shock, an d it flows around the cloud. In our interpretation 



of the velocity map (ISoker fc Kashill2012bl ) the post-shock gas flowing away from (towards) 
us forms the red-shifted (blue-shifted) component. Because of the inclination by 42° of the 
equatorial plane (the line from the center of the grid to the large WBE cloud in our model) 
to our line of sight the blue and red-shifted components are not symmetric, as we show in 
the next section. 

For comparison with observations later on (Section [5]) we build another very simple 
model for the structure of the WBE. In the second model for the WBE a smaller cloud is 
added to the cloud in the one-cloud model, as shown schematically in Figure [H The smaller 
cloud is denser than the large cloud, and is located closer to the primary star (for more 
details see Section |3]). The resulting density and velocity map for the two-clouds model is 
shown in Figure |H 

The pressure and temperature maps of both models are shown in Figure O It is 
evident from the temp erature maps that the post-shock primary wind cools very rapidly 



flSoker Sz Kashill2012bl ). As the iron responsible for the [Fe III] line is ionized, only the warm 
regions (T > 7500 K) contribute to the emission. However, it should be kept in mind that 
we do not include any radiative heating or ionization by the stellar radiation, and so we 
cannot quantify the contributions of the different warm regions. The pressure map shows 
the high pressure of the clouds that divert the primary wind. 



5. COMPARISON WITH OBSERVATIONS 



In comparing the velocity maps of our models to observations we are not only limited 
by the lack of knowledge of the mass distribution in the WBE, but also by the fact that 
the ionizatio n stage of iron and the ionizing radiation from the secondary star are poorly 



determined ( Madura et al. 2012[ Soker fc Kashi 2012b). T he variation of the velocity maps 



with orbital phase (IMadura et al.l l2012l : 



Gull et al 



201ll ) testify to the complexity of the 



problem. As our goal is to demonstrate the general velocity maps expected from the inter- 
action of the primary wind with the WBE, we use a very simple prescription. We assume 
that only gas with temperatures of T > 7500 K contributes to the [Fe III] line and that the 
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Fig. 4. — Same as Figure [3], but for the two-clouds model. The smaller cloud is in the bottom 
left part of the large cloud. 



-VI- 





Fig. 5. — The pressure (left column) and temperature (right column) maps for the one-cloud 
(top row) and two-clouds (bottom row) models. The scale in the pressure maps is logarithmic 
in units of dyn cm~^. The temperature maps are in K. 
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the emissivity goes as the density square (p^). This very simple approach should be kept in 
mind when examining our comparison with observations. For example, as we do not consider 
the ionization radiation from the secondary star, we have no dependence on orbital phase. 
Because of these uncertainties in our modeling, we basically make a qualitative comparison. 

To obtain the velocity components distribution maps, we integrate along the line o f 



sight for three velocity bins. The velocity bins are based on the bins in iGuU et al.l fl201l[ ) 
(-400 to -200 km s"\ -90 to 30 km s"^ and +100 to +200 km s"^), but we extend the 
boundaries of the three velocity bins by 30 km s~^ to each side, to account for the random 
velocity of gas in the post-shock regions. 



We compare our results with the observations of iGuU et al.l (120 llh in Figure [61 The 



ob servations at thr ee orbital phases are shown in the first three rows (taken from Figure 2 



of IGuU et al.ll201ll ). Orbital phase 12.0 is the spec troscopic event of January 2009, and the 



orbital period is 5.53 years ( iDamineli et al.l 120081 ). The fourth and fifth rows of Figure [HI 



show the synthetic velocity maps of our one-cloud and two-clouds models, respectively. The 
intensity in the observed maps (top three rows) is in square root scale, while the intensity 
in the synthetic maps (bottom two rows) is in logarithmic scale. As we do not follow 
the ionization degree of iron and the ionizing radiation, we cannot compare the absolute 
intensity values in our simulations to the observed values. Hence, we limit the discussion to 
a qualitative comparison between the spatial morphologies of the observed and the simulated 
maps. 

The overall morphology, both in the observed and simulated blue-shifted components 
(left column) is along the NE-SW direction, crossing the center but with some extension to 
the SE (lower left in the figure). The zero- velocity component (middle column) results mostly 
from the warm (T > 7500 K) WBE and slow wind around it, and reproduces quite well the 
qualitative structure of the observed velocity map. We see that even our simplified choice 
of initially spherical WBE in the one-cloud model manages to qualitatively reproduce the 
observed zero-velocity map. The red-shifted component (right column) in our simulations 
recreates the asymmetrical red-shifted component that appears in the NW region of the 
observed velocity map. The qualitative agreement with observations of the red component 
is better in the two-clouds model, as can be seen in the bottom-right panel of Figure |6l The 
small differences between the one-cloud and two-clouds models demonstrate that by different 
reconstructions of the WBE we can change the synthetic velocity maps. In principle we can 
find very good matching with observations. This is beyond the goals of the present paper 
as the parameters space of the WBE modeling is huge. Overall, we consider the qualitative 
agreement to be satisfactory considering the mentioned uncertainties and unknowns of our 
modeling procedure. We also compare the results of our simulations with the HST/STIS 
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Fig. 6. — .Top 3 rows show the observed velocity components at three orbital phases 
(iGuU et al.ll201ll ). with intensity in square root scale. The two bottom rows are synthetic 
maps of our simulations, with intensity in logarithmic scale. The fourth row from the top 
shows the velocity maps obtained for the single-cloud WBE model. The bottom row shows 
the velocity maps for the two-clouds WBE model (see Section [3]). The velocity bins in the 
simulated maps are extended to each side by 30 km s~^ to take into account random motion 
of the post-shock gas. 
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spatially resolved [Fe III] spectra of rj Car, as given by iMadura et al.l ((2012) • We reproduce 
the spectro- images for various position angles (PA) of the HST/STIS relative to the direction 
of the North in i] Car observations. The comparison is presented in Figures [7] and [8], and, as 
before, is only qualitative in nature. 

The agreement between the synthetic results and the observed line profiles for PA=+37.7° 
(see Figure [7]) is poor. However, considering the unknowns discussed above we cannot ex- 
pect a good match. None the less, we can identify observed features in our synthetic maps, 
as marked by letters in Figures [7] and [HI The blue-shifted velocities from km s~^ to 
—400 km s~^ are marked as feature B in Figure [71 The low- velocity components are marked 
as features A and C in Figure [7|). For PA = +68.7° (see Figure [S]), we reproduce the low 
velocity component present in the observed line profile at —40 km s~^ (marked as feature 
D in Figure [8]). We recall that in our simulations the velocity of the simple WBE model 
is km s~^, corresponding to the velocity value of feature D in the synthetic line profile. 
A more detailed and realistic construction of the WBE would give the WBE a velocity of 
~ —40 km s~^. There is some agreement between the general features of the observed and 
the synthetic line profiles for blue-shifted velocities (marked as features E and F in Figure 



Velocity [ km/s ] 



Velocity [ km/s ] 



Velocity (km/s) 



Fig. 7. — A qualitative comparison of the velocity along position angle PA = +37.7°. The 



right panel is the observed velocity profile in the [Fe III] line taken from (IMadura et al. 



20121 ). The color code in the synthetic maps is of the intensity calculated by integrating over 



along line of sight. Units are logarithmic and in relative values. Our synthetic profiles are 
given for the one-cloud and two-clouds models in the left and middle panels, respectively. As 
we do not calculate the ionization degree of iron and the ionizing radiation, the comparison 
is qualitative and limited to identifying large scale features. Features A and C mark low 
velocity components. Feature B marks the blue-shifted velocity component. The middle 
panel shows the same line profile for the two-clouds model. Although the two models seem 
to be almost similar, the small differences demonstrate that with correct construction of the 
WBE the velocity profile can be changed to better match observations. The emission from 
the center of the binary system (within 0."05) is not modelled in our simulations, and is 
covered by a white line. 



PA =68.7 , one-cloud model 



PA =68.7 , two-clouds model 
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Fig. 8. — Same as Figure [71 but for PA = -|-68.7°. Feature D marks the low velocity 
component. Features E and F mark blue-shifted velocity components. 
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SUMMARY 



We perform 3D gas flow simulations of the interaction of the r] Car primary wind with 
the dense gas clouds (WBE) in the equatorial plane of the binary system. Considering the 
complex and vaguely resolved structure of the WBE (see Figure [1]) , we choose two simple 
models for the WBE: (1) The WBE is modelled by a single spherical cloud; (2) The WBE 
is modelled by two spherical clouds, one la rge and one sma ll (see Figure [2] for more details). 
For comparison with [Fe III] observations (IGuU et al.ll201ll ) we construct synthetic emission 
maps at different velocity intervals (see Figure E]) • The exact reconstruction of the synthetic 
maps is made difficult by the high uncertainties in the ionization degree of iron in the WBE, 
and in the varying amount of ionizing radiation produced by the secondary star. Taking 
these into account we use a simple prescription to calculate the emission. We assume that 
the emission is produced only by warm gas with temperatures of T > 7500 K, and take 
the emission intensity to be proportional to p^. This set of simplifying assumptions allows 
only qualitative comparisons. As can be seen in Figure El we can qualitatively explain the 
velocity maps in the red, zero, and blue-shifted velocity intervals. The observed zero- velocity 
component is explained by the presence of slow dense warm material that constitutes the 
shocked wind and the WBE. 

We also compare the results of our simulation s with the HST/STIS spatially resolved 
[Fe III] spectra of 77 Car, as given by Madura et al.l (120121 ). While we are able to reproduce 
some of the large scale features in the observed line profiles, the overall agreement between 
the observed and the synthetic line profiles is somewhat limited. To obtain better agreement 
the WBE structure must be better resolved, and the ionization degree of iron and the ionizing 
radiation must be calculated. These are complicated tasks that are much beyond the present 
study. 

Based on the results of our simulations we conclude that the observed emission maps 
and line profiles can be explained as a result from the interaction betwee n the primary star 
wind and the WBE, as previously suggested by fISoker fc Kashil l2012b[ ). As the WBE is 
closer to the observer than the center of the binary system is, it follows from our results 
that the primary star (which blows the freely expanding wind) is also closer to the observer 
than the secondary star for the majority of the binary period, e.g. the longitude angle of the 
binary system is a; ~ 90°. Our results add to a growi ng number of obs ervat ions in favour 



of this longitude angle, and dispute the claim made by iGull et al.l ( 1201 ll ) and Madura et al. 



(I2OI2I ) that the observed velocity maps and line profiles exclude a value of a; ~ 90°. Our 
result shows that this orientation can account for the observations, and that the collision 
of the primary wind with the WBE cannot be ignored. We admit though that the winds 
collision should also be considered. Including both winds collision and the WBE interaction 
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is very demanding numerically, and should be the goal of future studies. 

Knowledge of the exact orientation of the rj Car binary system is of importance for 
correct interpretation of observations in the upcoming 2014 spectroscopic event. Knowing 
the correct value of lo will allow us to better understand the interaction of the binary system 
near periastron passages, and aid in further study of the accretion processes that took place 
during the Great and Lesser Eruptions, deepening our understanding of interacting LBV 
binary systems. 

We thank Amit Kashi for useful comments. This research was supported by the Asher 
Fund for Space Research at the Technion, and the USA- Israel Binational Science Foundation. 
The FLASH code used in this work is developed in part by the US Department of Energy 
under Grant No. B523820 to the Center for Astrophysical Thermonuclear Flashes at the 
University of Chicago. The simulations were performed on the TAMNUN HFC cluster at 
the Technion. 
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